Semi-analytical approach for the Vaidya metric in double-null coordinates 
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We reexamine here a problem considered in detail before by Waugh and Lake: the solution of 
spherically symmetric Einstein's equations with a radial flow of unpolarized radiation (the Vaidya 
metric) in double-null coordinates. This problem is known to be not analytically solvable, the only 
known explicit solutions correspond to the constant mass case (Schwarzschild solution in Kruskal- 
Szekeres form) and the linear and exponential mass functions originally discovered by Waugh and 
Lake. We present here a semi-analytical approach that can be used to discuss some qualitative 
and quantitative aspects of the Vaidya metric in double-null coordinates for generic mass functions. 
We present also a new analytical solution corresponding to (l/t;)-mass function and discuss some 
physical examples. 



PACS numbers: 04.40.Nr, 04.20.Jb, 04.70.Bw 

I. INTRODUCTION 

The Vaidya metric'l'l, which in radiation coordinates 
{w, r, 9, (j)) has the form 



ds^ 



2m(w)\ o T „T , , 

1 ^—!-]dw^ + 2cdrdw-^r^dn^, (1) 



where dQ"^ = d6'^+ sin 



c = ±1, is a solution of Ein- 



stein's equations with spherical symmetry in the cikonal 
approximation to a radial flow of unpolarized radiation. 
For the case of an ingoing radial flow, c = 1 and m{w) 
is a monotone increasing mass function in the advanced 
time w, while c — —1 corresponds to an outgoing radial 
flow, with m{w) being in this case a monotone decreasing 
mass function in the retarded time w. For many years, 
it has been used in the analysis of spherical collapse and 
the formation of naked singularities (For references, see 
the extensive list of 2] and also ^). It is known that 
Vaidya metric can be obtained, by taking appropriate 
limits in the self-similar case, from the Tolman metric 
representing spherically symmetric dust distribution _4J. 
This result has shed some light on the nature of the so- 
called shell- focusing singularities^, discussed in details 
in [1, la 13 ■ The Vaidya metric has also proved to be 
useful in the study of Hawking radiation and the process 
of black-hole evaporationj^, [ifl Ull H^ ll^. and more 
recently, in the stochastic gravity program[ij . 

Motivated by the longly known fact that the radia- 
tion coordinates are defective at the horizon[l3, implying 
that the Vaidya metric (Q) is not geodesically complete 
(see |la] for considerations about possible analytical ex- 
tensions), Waugh and Lakejl^l considered the problem 
of casting the Vaidya metric in double-null coordinates 



(m, V, 9, (f)). As all previous attempts to construct a gen- 
eral transformation from radiation to double-null coordi- 
nates have failed, they followed Synge[i3| and considered 
Einstein's equation with spherical symmetry in double- 
null coordinates ab initio. The spherically symmetric line 
element in double-null coordinates is 



rfs^ = -2f{u,v)dudv-\-r^{u,v)dn'^, 



(2) 



where f{u,v) and r{u,v) are non vanishing functions. 
The energy-momentum tensor of a unidirectional radial 
flow of unpolarized radiation in the eikonal approxima- 
tion is given by 



T„ 



1 



h{u,v)kah, 



(3) 



where ka is a radial null vector. As Waugh and Lake, 
we will consider, without loss of generality, the case of 
the flow along w-direction: ka = (0,1,0,0). The case 
of simultaneous ingoing and outgoing flows was consid- 
ered in [13. Einstein's equations are less constrained 
in such case, allowing the construction of some simi- 
larity solutions. In our case, the metric |(5J with the 
energy-momentum tensor Q reduce to the following set 
of equations [13: 



f{u,v) = 2B{v)dur{u,v), 

2A{v) 



dyr{u,v) 
h{u, v) 



-B{v) 1 



r(u, w) 



r^{u,v) 



(4) 
(5) 

(6) 



where B(v) and A{v) are arbitrary functions obeying, 
according to the weak energy condition. 



B{v)A'{v) < 0. 



(7) 



'Electronic address: lasaaQime. unicamp . br| 



Note that from Q), the regularity of f{u,v) requires 
B{v) y^ 0. Schwarzschild solution in the Kruskal-Szekeres 



form corresponds to ^' = 0. For A' ^ 0, the choice 

A' 



2B{v) 



\A'\ 



(8) 



allows one to interpret A{v) as the mass of the solution 
and V as the proper time as measured in the rest frame at 
infinity for the asymptotically flat case[l3- We assume 
here A{v) > 0. The radial flow is ingoing if A' > 0, and 
outgoing HA' < 0. Note that if the weak energy condi- 
tion (0 holds, the function A{v) is monotone, implying 
that the radial flow must be ingoing or outgoing for all 
V. It is not possible, for instance, to have "oscillating" 
mass functions A(v) in double null coordinates. 

The problem may be stated in the following way: given 
the functions A{v) and B{v), Eq. Q shall be inte- 
grated and r(u,v) is obtained. Then, f{u,v) and h{u,v) 
are calculated from I^J and ©. The arbitrary function 
of u appearing in the integration of (0) must be cho- 
sen properly to have a non-vanishing /(w, v) function 
from Q). Unfortunately, this procedure is not analyt- 
ically solvable in general. Waugh and Lake, neverthe- 
less, were able to find regular solutions for Eqs. iQ- 
(|H|) for linear {A{v) = \cv) and a certain exponential 
{A(v) = 4 (aexp(/3cu/2) + 1)) mass functions (A, a, and 
f3 are positive constants, c — ±1, corresponding to in- 
going/outgoing flow, respectively.). To the best of our 
knowledge, these are the only varying mass analytical 
solutions obtained in double-null coordinates so far. We 
notice, however, that Kuroda was able to construct a 
transformation from radiation to double-null coordinates 
for some particular mass functions 10, llj. 

In the following section, we will present a semi- 
analytical approach to attack the problem of solving Eqs. 
Q-® for general mass functions obeying the weak en- 
ergy condition 0. The approach allows us to construct 
qualitatively conformal diagrams, identifying horizons 
and singularities, and also to evaluate specific geomet- 
ric quantities. Before, however, we notice that one can 
solve analytically Eqs. Q-® also for the case of 



A{v) = -, 



(9) 



being k a massive parameter, and B — 1/2. With the 
mass function lO, Eq. ((SJ can be integrated as 

2 /"^ 2 

L(r,w,u) = P(M)w + 2Ke« -V2v e' ds ^ 0, (10) 

Jo 



where 



2\/2Vk k 



(11) 



and P{u) is an arbitrary function. It is quite easy to 
show that L{r, u,v) = is an invariant surface of (O, i.e, 
one has 



dL dL^ dL 1 

dv or ov V 



(12) 



along the solutions of (0 . From Q , we have 
/ = -kP'{u)v- 



2r 



(13) 



The function P{u) must the chosen to preclude zeros of 
f{u, v). With P{u) = —u/k, we have from (|10|l 



=2e^ 



r2^-t 

K Jo 



ds. 



(14) 



Equation (|14|l defines r{u,v). Its graphics and the corre- 
sponding conformal diagram is presented in Fig. ^ 
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FIG. 1: Constant r/n curves for the case A(v) = 
K^/v. The values of r/w shown are, from bottom to top, 
0, 0.5, 0.6, 1, 1.5, 2, 2.5, and 3. Future is to right and up. The 
corresponding conformal diagram is inserted. Here, r = is 
a naked singularity. 

The only spacetime singularity for t; > is r = 0. We 
recall that he Kretschmann scalar for metrics of the form 
(0) obeying 0-® is given hv\ii\ 






(15) 



The resulting spacetime corresponds to a naked singu- 
larity that vanishes smoothly, giving rise to an empty 
spacetime. Its causal structure is identical to the first 
linear case (A > 1/16) of Waugh and Lake, with time 
reversed (See Fig. |31 below). 

We finish this section noticing that algebraic manip- 
ulation programs (Maple, for instance) are able to find 
solutions of for the case A{v) oc w^ in terms of Airy 
functions, but they seem too cumbersome to be useful. 
Furthermore, their causal structure is always similar to 
the (A > 1/16) linear case of Waugh and Lake. Also, 
Kuroda has considered before the case of M{w) = fiw" 
for small w and n > 1 in radiation coordinates, obtaining 
some local properties of the singularity r = w = 0|lll]. 



II. SEMI-ANALYTICAL APPROACH 

Eq. Q along u constant is a first order ordinary dif- 
ferential equation in v. One can evaluate the function 
r(u, v) in any point by solving the u-initial value problem 
knowing r{u, 0) . The trivial example of Minkowski space- 
time (^ = 0), for instance, can be obtained by choosing 
r(u,0) = u/2. 

The curve r = 2A{v) is the frontier of two regions of 
the (w, r) plane where the solutions of 10) have qualitative 
distinct behaviors. Suppose, for instance, A'{v) > (and 
B = —1/2. The case of an outgoing radiation flow follows 
in a straightforward manner). For all points of the plane 
{v,r) below this curve, r' < (See Fig. |2J). Hence, any 
solution entering in this region will, unavoidably, reach 
the singularity at r = 0, with finite v. Suppose a given 
solution ri{v) with initial condition ri(0) = n enters into 
this region. As the uniqueness of solutions for ^ is guar- 
anteed for any point with r 7^ 0, solutions never cross 
each other in the plane {v, r) for r 7^ 0. Hence, any solu- 
tion starting at r(0) < n is confined the the region below 
ri(v), and it will also reach the singularity at r = with 
finite v. On the other hand, suppose that a given solu- 
tion rc{v) with initial condition re(0) = r^ never enters 
into region bellow the curve r = 2A{v). Any solution 




FIG. 2: In the region below the generic monotonic function 
r = 2A{v) (the apparent horizon), all solutions have r' < 0. 
Any solution that enters into this region will reach the singu- 
larity at r = with finite v. In the other hand, solutions con- 
fined to the r' > region always escape from the singularity 
and reach X^ . In this case, supposing that the mass function 
A{v) has an asymptotic value such that ri{v) < 2A{v) < ri{v), 
there exists an event horizon (the dashed line on the con- 
formal diagram) somewhere between the solutions ri{v) and 
ra{v). The vector T = dv ~ du points to the future for the 
ingoing radiation case. 



causal structure of the spacetime from their qualitative 
behavior and an appropriate choice for the initial condi- 
tions r{u, 0). For instance, in the case of Fig. [3 suppos- 
ing that the mass function A(v) has an asymptotic value 
such that ri{v) < 2A{v) < ri{v), an event horizon shall 
be located somewhere between ri{v) and rc{v). 

Event horizons are global structures. One needs to 
know the total (f — > (X)) spacetime evolution to define 
them. On the other hand, apparent horizons are defined 
locally!'2C|. Indeed, they are the real relevant ones in our 
semi- analytical approach. The curve r = 2A(v) is an 
apparent horizon. For a given w, all solutions in Fig. [21 
such that r{v) < 2A{v) will be captured by the singu- 
larity. Solutions for which r{v) > 2A{v) are temporally 
free, but they may find themselves trapped later if A{v) 
increases. The event horizon is the last of these solutions 
trapped by the singularity. 

We remind that Eq. Q requires that dur{u, v) 7^ 0, 
implying that, for v constant, r(u+du, v) ^ t(u, v) for all 
u. This conditions holds everywhere we have uniqueness 
of the solutions of (0 , provided that the initial conditions 
are such that 9„r(w, 0) 7^ 0. Moreover, if 9„r(M, 0) > 0, 
then dur{u, v) > 0, implying from Q that the sign of 
f{u, v) will be determined by B{v). For the ingoing radi- 
ation case, one has B{v) = —1/2, and from jSJ one sees 
that the vector field T = 3„ — 9„ is timelike. Our con- 
formal diagram are oriented such that this vector field 
points upward. If the condition dur{u, 0) 7^ holds, the 
conformal diagrams obtained for a given A{v) from differ- 
ent initial conditions are always equivalent. Some explicit 
examples will clarify the proposed approach. 



A. Linear mass function 

The linear mass function solution discovered by Waugh 
and Lake corresponds to (the ingoing radiation case, the 
outgoing one also follows in a straightforward way) A = 
Xv and B — —1/2. Note that with these choices, Eq. (Q 
does not satisfy the Lipschitz condition in r = 0, implying 
that uniqueness is not guaranteed for solutions passing 
there. This will be a crucial point to clarify the nature of 
the three qualitative different cases identified by Waugh 
and Lake, corresponding to A > 1/16, A — 1/16, and 
A < 1/16. We notice that the linear mass case was also 
considered in |^ in great detail and in a more general 
situation (the case of a charged radial null fluid) . 

The frontier of the region where all the solutions reach 
the singularity at r = is, in this case, the straight line 
r = 2Ar. Taking the w-derivative of Q, one gets 



1- 



2Xv 



- 1 



(16) 



starting at r(0) > r^ is, therefore, confined to the region 
above r^iv) and will escape from the singularity at ?' = 0. 
Keeping in mind that these solutions correspond to null 
trajectories with u constant, it is possible to infer the 



The regions in the plane {v,r) where the solutions obeys 
r" = are the straight lines 



^(i±VT^T6a) 



(17) 



One can reproduce the analysis of Waugh of Lake for the 
three quahtative different cases according to the value of 
A by considering the possible solutions of 1)17(1 . For this 
purpose, we will consider the solutions of with the 
initial condition r{u, 0) ex u. 

For A > 1/16, the case showed in Fig. ©, r" < 
for all points with v > 0, and the only relevant fron- 
tier is the r — 2Xv straight line (the apparent horizon). 
All solutions of ^ are concave functions and cross the 
r' = line, reaching the singularity with a finite v which 
increases monotonically with u. The causal structure of 




FIG. 3: For A > 1/16, all solutions are concave for v > 
0. Hence, all null trajectories along v necessarily reach the 
singularity at r = 0. The figure corresponds to the case A — 
1/10. The conformal diagram is inserted. 

the corresponding spacetime is very simple. There is no 
horizon, and the future cone of all points ends in the 
singularity at r = 0. 

Each curve in Fig O corresponds to a constant u slice 
of the the full solution r(u, v). It is possible, in principle, 
to reconstruct the 2-diniensional surface r(u, v) and to 
plot its lines of constant r on the (m, v) plane, as done 
in the Fig^ Q , reproducing all original results of Waugh 
and LakeflTj. However, these lines are not necessary to 
construct the conformal diagrams. 

For A = 1/16, r" ~ on the line r = v/A. This straight 
line itself is a solution. All other solutions are concave 
functions. We have two distinct qualitative behavior for 
the null trajectories along u constant. All solutions start- 
ing at r(0) > are confined to the region above r = w/4. 
They never reach the singularity, all trajectories reach 
I+. However, in the region below r = v/A, we have in- 
finitely many concave trajectories starting and ending in 
the (shell-focusing) singularity. They start at r(0) = 0, 
increase in the region between r = v/A and r — v/8, 
cross the last line and reach unavoidably r — again, 
with finite v. The trajectory r = v/A plays the role of an 
event horizon, separating two regions with distinct qual- 
itative behavior: one where constant u null trajectories 



reach X"*" and another where they start and end in the 
singularity. This situation is showed in the Fig. 0] This 
behavior is only possible, of course, because the solutions 
of (jSJ fail to be unique at r = 0. 




FIG. 4: For A = 1/16, we have two relevant lines in the (r, v) 
plane: r = v/A and r — v/8. The first one plays the role of 
an event horizon, separating the two regions of distinct qual- 
itative behavior for the null trajectories along v. Beyond the 
horizon, all solutions escape to infinity. Inside, all solutions 
reach the singularity r = with finite v. The conformal dia- 
gram is inserted. The singularity at r = is a shell focusing 
one. 

For A < 1/16 (Fig. 0), we have three distinct regions 
according to the concavity of the solutions. They are lim- 
ited by the two straight lines (|17() , that are also solutions 
of ^. Between them, solutions are convex. Above r+ 
and below r_, solutions are concave. The last line is the 
inner horizon. Inside, the null trajectories with constant 
u start and end in the singularity r = 0. The line r+ 
is an outer horizon. Between them, the constant u null 
trajectories starts in the naked singularity and reach X+. 
Beyond the outer horizon, all constant u null trajectories 
escape the singularity. 



B. Collapse of a radial null fluid 

As an example of the proposed approach applied to 
new mass functions, let us consider first the case of a 
collapse of incoming radial null fluid, starting with the 
empty space. This could be described, for instance, by 
the mass function 



Aiv) 



-(1 +tanhpv). 



(18) 



The corresponding spacetime is empty for v —f — cx3, it 
receives smoothly a radial flow of null fluid, and finishes 
as a black hole of mass m for w — > oo. The solutions of 
^ are presented in the Fig. With the mass function 
(|18|l . we always have a singularity at r = for f ^ — oo. 




FIG. 5: With A < 1/16, we have two horizons r_ and r+. 
Any solution starting above r+ always escapes to I+. Solu- 
tions starting inside the inner event horizon ri always reach 
the singularity at r = 0. Between the horizons, solutions 
starts in the naked singularity but escapes to X+. The figure 
corresponds to the case A = 1/20. The conformal diagram is 
inserted. 




FIG. 6: The case A{v) = m (1 + tanhp?;) /2. The solutions 
tend asymptotically to r=0 when v —> —oo. A zoom of the 
range — 3 < u < — 2 illustrating such asymptotic behavior is 
inserted. In this case, r = is always singular for v > — cxd. 
This plot corresponds to mp — 1, the situation is qualitatively 
identical for other cases. The conformal diagram is inserted. 



The constant u null trajectories approach r — asymp- 
totically for V — > — cx). Some of them, confined by the 
horizon, are captured by the collapse and reach the sin- 
gularity, and others scape to T+. Any initial condition 
obeying 9,jr(u,wo) 7^ leads to the same causal struc- 
ture. 

The proposed semi-analytical procedure can be used 
also to derive some qualitative results. One could, for 
instance, locate accurately the event horizon in this case. 



Since we do not expect deviations from the Minkowski 
spacetime for w — > — oo, we could start with initial con- 
ditions such that r(u, —cxd) = u/2. Then, we identity 
in the solutions of (jSJ that one (the horizon rh{v)) ap- 
proaching asymptotically the curve r — 2A{v) for large 
V. For large values of vq, the horizon would be located 
at (2rh(-^;o),-i;o). 



C. Black hole evaporation 

As a last example, let us consider the physically rel- 
evant case of black hole evaporation. We remind that 
much work has been devoted on Vaidya metric to describe 
black hole evaporation 9, 10> lllLll2l Il3|. Our approach 
allow us to analyze this question in the more convenient 
double-null coordinates. The process of black hole evap- 
orations supposes A' < and an outgoing radiation flow. 
Hence, one needs to chose B{v) = 1/2 in our procedure. 
With this choice, we see from Q that f{u, v) > 0, im- 
plying that the field dy — du is not timelike anymore. It 
is namely its orthogonal vector field d^, + du that will 
define the time evolution in this case. In order to keep 
the future direction pointing upward in the conformal di- 
agrams, one needs to re-orient the coordinates u and v. 
The appropriate transformation here is 

V ^ u, u —f —V. (19) 

With this choice and B(v) — 1/2, equation (O reads 



^ . ^ 1 / 2A( 
dur[u,v) = 2 ' 1 



r{u,v) 



(20) 



which is formally identical to the equation for the ingoing 
case characterized by the mass function A{—u) with radi- 
ation flow along the u direction. Hence, the description of 
an outgoing situation with decreasing mass function A{u) 
can be obtained by applying the transformation (|19|l for 
the ingoing case corresponding the mass function A{—v). 
With the help of ifT^ . one could construct the confor- 
mal diagram corresponding to the mass function A{u) = 
m{l — tanhpu)/2 starting from the last example. Such 
a case, however, does not correspond really to an evapo- 
rating black hole since, according to this mass function, 
the black hole never disappears, the singularity at r = 
is ever present for m < oo. In a real process of black 
hole evaporation due to Hawking radiation, the black 
hole loses mass with a rate M oc Af~^, disappearing 
completely in a finite time and giving rise to an empty 
spacetime. The mass function 



A{u) 



— TTitanhpu, u < 0, 
0, u>0. 



(21) 



can be used to simulate the vanishing of a black hole 
in a finite time. For the ingoing radiation case A{—v), 
the situation is similar to the last example, but with the 
crucial difference that for u < the spacetime is empty, 



and the solutions can reach (and cross) r = 0. We no- 
tice that the form of A{u) is not important. In order to 
construct the conformal diagram of a evaporating black 
hole, one only needs that A{u) — after some u < oo. 
In this case, the full conformal diagram (see Fig. [TJ is 




FIG. 7: Conformal diagram for a evaporating black hole with 
mass A{u) given by 12111 . The region I corresponds to the 
black hole, and can be obtained, by using the transformations 
H19|l . from the ingoing radiation case. Region II is the empty 
space left after the vanishing of the black hole at u = 0. The 
event horizon is the line EH. The dashed CH line is a Cauchy 
horizon. Beyond it, we have breakdown of predictability due 
to the naked singularity left by the black hole. 

constructed adding to that one obtained from the ingoing 
collapse case by using the transformation H19|l . a patch 
of Minkowski spacetime, just after the instant u = 



corresponding to the vanishing of the black hole. The 
resulting conformal diagram represents a black hole that 
evaporates leaving behind a naked singularity. 



III. CONCLUSION 

We presented a semi-analytical approach to construct 
the Vaidya metric in double-null coordinates for general 
monotonic mass functions. We applied it to elucidate the 
nature of the three distinct qualitative causal structures 
identified by Waugh and Lake|17j for linear mass func- 
tions, to construct a solution corresponding to a smooth 
radial collapse of a null fluid, starting from empty space 
and finishing with a black-hole, and to construct the con- 
formal diagram of a evaporating black hole. 

The approach involves an arbitrary function r{u,0). 
However, any choice for which dur{u, 0) 7^ will result 
in causally equivalent regular (/(m, v) 7^ 0) space-times. 
The analytical solutions of this problem also involve an 
arbitrary function, as P{u) in (|1U|) . that must be chosen 
properly in order to get a regular spacetime. 

The double-null form of the Vaidya metric is suitable 
to the study of quasi-normal modes of time depending 
solutions. In double-null coordinates, the equations for 
scalar perturbations are separable and we obtain an ef- 
fective Schrodinger equation with a time depending po- 
tential for the perturbations. This situation has been 
recently considered in the heuristic analysis of [2]|. The 
full problem is now under investigation. 
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